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Q_i' We explore a new class of axion models in which some, but not all, of the 

left-handed quarks have a Peccei-Quinn symmetry. These models are poten- 
tially afflicted by flavour changing neutral currents. We derive the bounds on 



X 

the Peccei-Quinn symmetry-breaking scale from bounds on the K + — > 7r + a 
branching ratio, showing that in some cases they are even stronger than the 
astrophysical ones, but still not strong enough to kill off the models. 
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I. INTRODUCTION 



One of the most persistent problems in particle physics is the so called strong CP problem. 
The CP invariance of the strong interactions can, in principle, be spoiled by the addition of 
the allowed '0-term' 

However, experimental limits on the neutron dipole moment imply that 9 < 10~ 9 PJ3L This 
unnatural smallness of the ^-parameter is precisely the strong CP problem. 

Over the years many solutions have been proposed for resolving this puzzle. A most 
appealing idea is the one brought forward by Peccei & Quinn 0]. They postulated the 
existence of an extra global anomalous U(l) symmetry. The accomodation of the new 
charges requires, at least, one extra Higgs doublet. With the help of the PQ-symmetry 9 
becomes dynamical and is driven to zero. However, since the new symmetry is not manifest 
in our world, it has to be spontaneously broken. As a result of the Goldstone theorem, a 
pseudoscalar particle appears [|], called the axion. Although one would expect the axion 
to be massless, since it is the Nambu-Goldstone (NG) boson of a global symmetry, it is 
actually not. The reason is that, as mentioned before, the PQ-symmetry is an anomalous 
one, spoiled by instanton effects, a fact that forces the axion to pick up a small mass through 
the axion-gluon-gluon anomaly. 

In principle, there is a plethora of axion models as there are many ways of assigning the 
PQ charges to the quark fields, and a lot of freedom to introduce extra Higgs fields. The 
original model ||, where all quarks were assigned the same charge and where two Higgs 
doublets were used, was experimentally ruled out. In order for the Peccei-Quinn solution 
to survive, another class of axion models was invented, called the invisible axion models, 
where a Higgs singlet was added H[j|. As a result, the axion decay constant v a became 
much higher than the electroweak scale and this enabled the axion coupling to matter to 
become much smaller than those of the original model. However, the value of v a cannot 
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be arbitrary. It is constrained from below by astrophysical observations and from above by 
cosmological considerations. The current values on these limits are 10 10 < v a < 10 11_12 GeV 
leaving a small window for the axion to exist. The value of the upper limit depends 



on the cosmological scenario one favours, that is either inflation or cosmic strings, the latter 
being the most restrictive. 

The two mostly talked about axion models are the KSVZ || and the DFSZ [frj. This does 
not mean that they are the only ones allowed, as we have demonstrated in a recent paper 



in which we explored the consequences of assigning different PQ charges to different 



right-handed quarks. In this paper we allow the left-handed quark doublets to have different 
PQ charges, which forces us to confront the problem of flavour changing neutral currents 
(FCNCs). Drawing on the work of Feng et al fllf , we find that FCNCs constrain the axion 
scale even more strongly than the astrophysical arguments (with one exception). The most 



stringent limit comes from the decay K + — > ir + a ||15|| . There are of course other bounds on 
the axion scale, for example /i — > ea [[□]]. However, these turn out not to be as strong as 
the K + — > ix + a constraints, so we do not consider them here. 

The class of models we consider represent a minimal digression from the DFSZ model, 
in that instead of all the left-handed quarks having the same PQ charge, only two of them 
do. That requires further extension to the Higgs sector, by adding at least one (two in the 
super symmetric case) new Higgs doublet (s). As a consequence there are up to four other 
neutral scalars and up to three other pseudoscalars in our theories, which could contribute 
to CP-violating processes, in particular neutral meson mixing (a fifth scalar coming from 
the singlet is irrelevant for low energy applications). Limits on the mixing translate to limits 
on the masses of these pseudoscalars of about a few hundred GeV to a few TeV, which in 
turn bound parameters in a complicated three or four-Higgs potential. Such potentials are 
somewhat problematic in general, as there is a severe hierarchy problem to solve in order to 
separate the electroweak and Peccei-Quinn breaking scales. 
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II. DESCRIPTION OF THE MODELS 



The main idea in this paper is to assign appropriate PQ charges to the Higgs fields and 
consequently to quarks, so that FCNCs can be present at the tree level of the axion-quark 
couplings. Then, it will be possible to derive constraints on the axion decay constant based 
on recent experimental data from the absence of such processes. In order to take advantage 
of FCNCs, we give different charges to the left-handed sectors of some of the quarks also 
(unlike the DFSZ where all left and all right handed quarks have the same PQ charges). 
This is a special case of a minimal change from the DFSZ. The general consistency rules 
that govern such digressions are the topic of future work |H| . For constructing these models 



one needs at least three doublets n and one singlet 0. If we allow four Higgs doublets, we 
have the possibility of making the model supersymmetric, as discussed below. The general 
structure of the Yukawa couplings is 

£Y = f% U {qLi<Pn u U Rj ) + f% d {qLi<Pn d d' Rj ) + h.C (2) 

where n u — 1, 3, n& = 2, 4 and i,j = 1,2,3 are flavour indices. This will result in six different 
models depending on which quarks have PQ charges. Following the notation of our previous 
paper fll2"| , in the first three the 'special' doublets are either the (u,d) L , or (c,s)x, or (t,b)i,, 
labeled by I, IV, II, and in the last three ones either {u,6)l and (c,s)x, or (u,d)i and (t,b)i, 
or (c,s)l and (t,b)i respectively, labeled by V, III, VI. In the absence of supersymmetry, 
where the appearence of 03 is forbidden as the superpotential must be holomorphic, one can 
put 04 = 03 = 2^203. Note that all quarks are left as flavour eigenstates for the time being. 
The general lagrangian, part of which is the Yukawa sector mentioned above, possesses a 
PQ symmetry. The most general PQ transformations are 

u'm — > e iaT ^u' Rj 
d'm — e^d' Rj 

q'Li ^ ^q' Lj (3) 
n _ e <Q»«0 nj n= 1,2,3,4. 



where Q% = Qi = 1 and Q3 = Qa = 0- The transformation matrices for the left-handed 
quarks, T L , and the right-handed u-type quarks, T Ru , are listed in Table | for every model. 
In the case of the right-handed d-type quarks the transformation matrices Tr^ = in all 
cases thus not listed in the Table. As an example, let us consider Model I, the tranformation 
matrices of which are 
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Rd 
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Their application fixes the Yukawa couplings f^' to have zeros in certain entries. For the 
■u-type quarks the relevant Yukawa matrices are /?• and /£. In the first one, the only non- 
zero element is and in the second we must have = /j 5 - = 0. On the other hand, for the 
d-type quarks, /?■ and being the appropriate Yukawa matrices, we need /f ■ = = 0, 



and f 4 



0. 



Our next step is to determine the axion decay constant in terms of the vacuum expecta- 
tion values of the Higgs fields and the mixing with the Z°. Suppose that a' and Z are the 
would-be Goldstone bosons before instantons are taken into account. We define a' to be the 
massless axion and Z the longitudinal degree of freedom of the Z° boson. Suppose also that 
a and az are the angles conjugate to the PQ and Z transformations, so that 
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exp [i(2a — az)], 



exp (—laz) 
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exp [i(a + a z )\-, 



exp \iaz) 



(4) 



72 Veiq> 2 a 

where v%, V2, V3, v±, v are the vacuum expectation values of the Higgs fields. In order to 
separate the axion from the Z°, one comes to the following equation 

( -\ 



(A 



V a 

Via v n 



\ az ) 



(5) 



where the 2x2 matrix on the right hand side of fl5|) is the most general one compatible 
with this requirement. A simple comparison of (||) with the kinetic terms of the NG bosons 
(coming from the kinetic terms of the Higgs fields) 0] gives the expression for the axion 
decay constant and for the electroweak breaking scale 



9 a (Av! + vl)(v! + vl) + 9vW 2 

''" = 4 V + v\ + v\ + 4 + vl ' {b) 

vew = v u = ^vl + vl + vj + vl = 246 GeV, (7) 

vio = 2 1 . (8) 

yjv% + v\ + vl + V% 

As we see, it has the same value for all six models. Furthermore, it is essentially equal to v, 
the PQ symmetry breaking scale. 



III. AXION-QUARK COUPLINGS AND INDUCED FCNCS 

Let us take now a closer look at the axion couplings to quarks. In the flavour basis the 
relevant term of the QCD lagrangian is 
d^nf 

Cim = -^K7,((l - l,)Ti + (1 + 75 )n J >; + <7,(1 - Ts)^^] ( 9 ) 
2v a 

It is possible to diagonalise, the generally non-diagonal, quark mass matrix by a bi-unitary 
transformation. In more detail, there are unitary transformations that relate the flavour 
basis with the mass one, of each of the u and d-type quarks 

u 'm = UrUrj, d Ri = D l ^d Rj} 
u 'u = U l ^u Lj , d' Li = D l [d Lj . 

Applying these transformations to (§) and going to the mass basis, the lagrangian takes the 
form 

£int = — [2ua^ 5 S l i u Uj - dap ( 1 ~ 7s) S'^dj] (10) 

where Sl u = U\T L U L and = D^ L T L D L . By definition, the Cabbibo-Kobayashi-Maskawa 
matrix is Vqkm = U^Dl, so it is obvious that 



Shd — V^ km Sl u Vckm, 



(11) 



thus being possible for FCNCs to be present in the d-type quark sector, both in the vector 
and in the axial- vector part of the Lagrangian. Furthermore, it is evident from the structure 
of the w-type Yukawa couplings, that Ul and Ur have a block diagonal form and thus 
Sl u = T L in all cases. So eq. ([O]) becomes 



SlcI — V ckm T l Vckm, 



Combining the data from Table | and eq. (|T2|) one finds 



(12) 







VuidVuiS 




T d = 




V* V 


v: iS v Uih 






V* V 





(13) 



for Models I, II, IV, where i = 1,2,3 labels the PQ-charged quarks and 



/ v; id v Uid + v* d v u . d v: id v UiS + v* d v UjS v: id v Uib + v* d v Ujb ^ 



\ 



(14) 



v: lS v Uld + v: ]S v U]d v*. s v UiS + v: ]S v U]S v: iS v Uib + v: jS v Ujb 

V u,b v u t d + v*. b v Ujd v*. h v UiS + v*. b v UjS v* ib v u . b + V*. b V Ujb j 

for Models III, V and VI, where i ^ j also labeling the relevant PQ-charged quarks. As we 
shall see in the following section the interesting part of the interaction lagrangian (|ToD is the 
one giving the transition of s — > d quarks. In this case 



int 



2v a 



[sin(g sd + 9sdl5)d+ h.c] 



(15) 



g^ d and gf d being, by definition, the vector and axial vector parts of the a — s — d coupling. 
Combining eqs. (0), (0), dH) and ( |i~5|) one finds 

Models I, II, & IV 

(16) 

V: id V UiS + V*. d V UjS Models III, V, & VI 

Concerning the axial coupling, as we shall see in the next section is of no importance, since 
only the vectorial one is involved in the calculation of the rate K + — > n + a. The values 



9sd 



V* V 



for the CKM elements used are [1T4| IK 



'ud\ 



0.98, \V US \ 



0.22, \V cd \ w 0.22, \V ct 



0.97, 



Wt 



td\ 



9.1 x 10~ 3 and \VJ w 3.9 x 10" 2 . 



IV. EXPERIMENTAL CONSTRAINTS 



As described in the previous section, the axion can take part in flavour- changing pro- 
cesses. One can extract lower bounds on the axion decay constant from experimental data 
concerning such processes. The tighter constraints come from transitions between the first 
two generations, whereas bounds involving the third one are much weaker. The processes 
that produce the most stringent limit are the ones involving rare K decays. The most 
suitable one for our discussion is K + — > ir + a, the decay rate of which is (TT1 

- n+a) = -^4^/(1 - 4) 3 |^i(0)| 2 (17) 
16% v l a m z K 

where -^1(0) is the form factor Fi(q 2 )(p + = (TC + (p')\s , y >J, d\K + (p)) at zero momentum 
transfer and is of order unity, being exactly 1 in the case of exact SU(3) flavour symmetry. 
The experimental data sets an upper limit on the branching ratio [|l5j Br(K + — > 7r + a) < 
3.0 x 10 -10 (at 90% confidence). This leads to a lower bound on the axion energy scale 

v a > 1.7 x 10 11 x gj s GeV (18) 

Taking into account eq. ( f[q) and the current values for the relevant CKM matrix elements 



14 1 , the above expression yields lower limits on v a . The results are listed in Table |TJ. 

It will be very instructive to compare these results with the astrophysical ones since the 
latter are so far considered to be the most severe. It is a well known fact that among the 
astrophysical limits the far more restrictive are the ones coming from SN1987A, bounding 
the axion- nucleon-nucleon coupling. The limit of this constraint, taking many body effects 



into account, is 16 



(/4 + 2/4) 1/2 <2.85xl0- 10 . (19) 



A similar analysis as the one performed in [12], normalising the PQ charges for N = 1 or 



N = 2 depending on the case, yields the following lower bounds on the axion decay constant 

v a > 0.35 x 10 10 x (A/i 2 -B/i+ C) 1/2 GeV (20) 



where 

_ Vw _ v% ~ v i 

\l — — 2 

2 2 

and — j-^- < /i < The values of the coefficients A, B and C are summarised in Table 

V EW V EW 

111] for each model. One can easily see from eq. (p0[) that the most stringent limits come 
for fx = —1, in the limit where V2 = v% = = 0. However, comparison of these results 
(plotted in Fig. [1]) to the ones coming from FCNCs reveals that for all models, except the 
third one, the latter constraints are more severe. Especially in the case of the fifth model the 
limit is almost up to 10 11 GeV, very close to the upper bound coming from the cosmological 
scenario of cosmic strings |I(J. The exception of the third model is due to the fact that 
VtdVts "C V* d V us ~ V* d V cs , also responsible for the similar behavior of the rest of them. 

V. CONSTRAINTS ON MASSIVE SCALARS AND PSEUDOSCALARS 

A theory with four Higgs doublets and a Higgs singlet has an additional five massive 
neutral scalars and three massive pseudoscalars. Apart from the scalar coming from the 
singlet, the rest can in principle mediate neutral meson mixing in the low energy sector. For 
example, BB mixing proceeds via the operators 

1 



O b BS = hUdb)- r (bd), 

OZ s = hUdjsb)j^(Hd), (21) 



- l" ;»»—, 

k 



for the scalar and the pseudoscalar case respectively, where h% d is the flavour- changing cou- 
pling and Mh k ,MA k are the masses of the k— th scalar and pseudoscalar. Evaluating the 



matrix elements for BB mixing in each case one finds [17| 



^ B m^fW B (n^^ + 1 ) (22, 

where and m# are the B decay constant and mass respectively and m^^vrid are the 
masses of the b and d quarks. These expressions were evaluated using the vacuum insertion 

9 



approximation in both cases, in order for a direct comparison to be possible. Of course one 
may argue that vacuum insertion, although a good approximation in the pseudoscalar case, 
is not as good in the scalar. For an order-of-magnitude estimate it should however suffice. 



In the heavy quark approximation, where mg ~ m& + n%d, this gives a mass splitting of [17 



2 f 2 
Am B ~ ^h 2 bd j^m B , 



h 

Am B ~ 4hl d ^m B , (23) 

for the scalar and the pseudoscalar interactions respectively, where we have dropped the 
k index, since we assume that the larger mass splitting is due to the lightest among the 
scalars and pseudoscalars. It is obvious that pseudoscalars provide stronger constraints by 
a factor of 6. In order to extract constraints for their masses we must first comment on the 



values hf,d can take. Cheng & Sher [JTSfl have argued that in a wide class of models where 



no fine-tuning is assumed, the flavour-changing couplings are of the order of the geometric 



mean of the Yukawa couplings of the generations involved, tif, = n^fifj- Thus in our case 



we have 

1,2 _ 2 f f _ 2 m b' m >d (OA\ 
n bd — K bdJbJd — K bd . ^Z4J 

{(f>n){<Pm) 

where n and m take the values 2 or 4 depending on the model. Since we do not know the 
expectation value of each Higgs field we will attempt an estimate of the order-of-magnitude 
for the Higgs masses. Taking the value for the decay constant of the B meson to be f B ~ 
170 MeV from lattice calculations ||19[ and the experimental value for Am B < 3 x 10~ 13 GeV 



20| , we obtain 

M h > k u x 1.3 TeV, M A > K hd x 3.2 TeV. (25) 

As we can see the results are in the TeV scale. This means that it is necessary only to 
make the not-so-unreasonable assumption that k's which involve the first generation can 
take values much less than one, in order that the constraints ( [25] ) can be further weakened 
and the masses can take values within the range of a few hundred GeV, a natural Higgs 

10 



mass for an electroweak theory. However, we are still left with a hierarchy problem, that is 
common to all invisible axion models; how to account for the wide range of scale between 
100 GeV and 10 12 GeV. Our model is no improvement is this regard. It should in principle 
be possible to make the models supersymmetric, when there are four Higgs doublets, in 
which case it is consistent to set the singlet-doublet couplings to be extremely small without 
fear of radiative corrections. 

VI. CONCLUSIONS 

It is evident from the above analysis that there is a lot of freedom in choosing the PQ 
charges in the quark sector. In this paper we studied a new class of axion models, where 
the left-hand sector of certain quark flavours (but not all) were assigned PQ charges. As 
a consequence, FCNCs are induced, which can be used to provide a lower bound on the 
axion decay constant. It was shown that for certain models the limits on these bounds 
are more severe than those coming from astrophysics, with the most striking example the 
case of Model V, although not severe enough to rule them out. We have also estimated 
constraints coming from neutral meson mixing induced by other scalars and pseudoscalars, 
which constrain their masses to be in the range of a few hundred GeV to a few TeV. 
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TABLES 



Model 


T L 




Tr u 




I 


diag (1,0,0) 




diag 


(-1,0,0) 


II 


diag (0, 0, 1) 




diag 


(0,0,-1) 


III 


diag (1,0,1) 




diag 


(-1,0,-1) 


IV 


diag (0,1,0) 




diag 


(0,-1,0) 


V 


diag (1,1,0) 




diag 


(-1,-1,0) 


VI 


diag (0,1,1) 




diag 


(0,-L-i) 


TABLE I. Transformation matrices for all left-handed and ti-type right-handed quarks, con- 


cerning 


the models discussed in the text. 


Tad = in all cases thus not listed. 


Model 


Charged 


Vector 




Axion scale 




doublets 


coupling 




(xlO 10 GeV) 


I 


ud 


v: d v us 




3.7 


II 


tb 






0.0061 


III 


ud, tb 


V ud V us + V t * d V ts 




3.7 


IV 


cs 


VcWcs 




3.6 


V 


cs, ud 


v: d Vus + v c * d v cs 




7.3 


VI 


cs, tb 


v: d v cs + v t * d v ts 




3.6 



TABLE II. Limits on the axion scale v a from the flavour-changing process K + — > n + a. 
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Model 


A 


B 


C 


I 


3.98 


2.47 


0.42 


II 


3.98 


1.37 


0.79 


III 


3.98 


4.01 


1.22 


IV 


3.98 


1.52 


0.61 


V 


3.98 


3.84 


1.31 


VI 


3.98 


2.29 


4.57 



TABLE III. Table of the values of the coefficients in Eq. (|20| ) for the axion models considered 



in the text. 
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